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■ Abstract. Irrotational dust spacetimes with vanishing magnetic Weyl curvature are 

O , called silent universes (Matarrese, Pantano and Saez, 1994). The silent universe 

conjecture (Sopuerta 1997, van Elst et al. 1997) states that the only algebraically 
general silent universes are the orthogonally spatially homogeneous Bianchi I models. 
\Q : In the same paper by Sopuerta this was confirmed for the subcase where the spacetimc 
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also admits a group G3 of isometries. However the proof contains a conceptual mistake. 
. We recover the result in a different way. 
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O ■ 1. Introduction 

O 

^\ Irrotational dust (u a = 0,p = 0) solutions of the Einstein perfect fluid field equations 
Sh . 

Rab - -Rgab = pU a U b (1) 

with vanishing magnetic part H ab of the Weyl tensor pQ were called silent universes for 
the first time by Matarrese, Pantano and Saez [2|. The reason for this nomenclatura 
was that signal exchange in GR can only occur via sound and gravitational waves, none 
of wich modes are allowed when p = H ab = 0. From a more mathematical point of 
view, the evolution equations for the matter density p, the three non-zero eigenvalues 
6*1? 6 I 2? ^3 of the expansion tensor dab an d the two independent eigenvalues (e.g. Ei,E 2 ) 
of the remaining electric part E ab = C acba u c u d of the Weyl tensor, form an autonomous 
system of ordinary differential equations (see also below); no spatial gradients appear, 
such that each fluid element evolves indeed as a separate or 'silent' universe, once the 
constraint equations are satisfied by the initial data. This setting looked very appealing 
towards numerical schemes and simulations in astrophysical and cosmo logical context, 
e.g. for the description of structure formation in the universe and the study of the 
gravitational instability mechanism in general relativity, where a clear motivation for 
taking p = u a = 0, H ab « was given in 

The only allowed Petrov types for silent universes are O, D and I: the Friedman- 
Robertson- Walker dust metrics exhaust the type class, while all Petrov type D 
solutions are known explicitly too. They are characterised jl] by the fact that the Weyl 
tensor is degenerate in the same plane as the shear and are given by the Szekeres [H] 
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family, including the subcase of the Ellis |B] LRS class II dust models and such well 
known examples as the Lemaitre-Tolman-Bondi model and the orthogonally spatially 
homogeneous Kantowski-Sachs model. For the algebraically general case (Petrov type 
I) the situation is rather different, the only known silent universe being the orthogonally 
spatially homogeneous Bianchi I model. 

In two independent papers by Sopuerta [7| and van Elst et al. the propagation 
of the constraint = was shown to give rise to a triplet of (in principle) infinite 
chains of equations, identically satisfied for Petrov type 0, D and spatially homogeneous 
silent models, but not in general. What happens is that for non-spatially homogeneous 
Petrov type I silent universes the initial values for the 6 a , E a and fi must lie on a (low- 
dimensional) algebraic variety, i.e., there exist severe polynomial relations between these 
variables. As it was conjectured by the authors of 0IE], there are good reasons to believe 
that the variety is actually empty, i.e., that there exist no non-spatially homogeneous 
Petrov type I silent models at all, at least for zero cosmological constant. This is the 
so called 'silent universe conjecture', no definitive proof of which has been given until 
now, mainly because the degree and the number of terms of the polynomials are massive. 
However, we mention here that in the case of strictly positive cosmological constant, the 
above variety has been shown to contain three (analogous) two-dimensional components, 
corresponding with two explicit new families of metrics each jHj. 

In the present paper, we correct a reasoning in [Ij and reconfirm the silent universe 
conjecture (for vanishing cosmological constant) in the subcase where the spacetime 
admits a group G3 of isometries. 

For silent universes, one can deduce in general (Proposition 3 of [Zj) that there 
always exists an orthonormal basis {e , e a } (with e = u) of common eigenvector fields 
of the expansion tensor a & and the remaining electric part = C ac bdU°u d of the 
Weyl tensor, which are all parallelly transported along the matter flow lines and are 
hypersurface orthogonal, such that a local coordinate system (t, ar ) exists for which 
eo = ^ and e a oc Below we will use the notation d a for e a whenever the frame 
vector fields act as differential operators. 

In jTj algebraically general silent universes were investigated in such coordinates, 
and in section 5 restrictions on the coordinate components of a generic Killing vector 
field K were obtained. However, in the discussion of Case (iii) a time dependence via 
K z in (108) was overlooked and therefore the conclusion that K l is constant and K x " 
only depends on x a does not necessarily hold. But even if this conclusion were valid, a 
coordinate transformation x a — > x a ' which transforms the components of a fixed Killing 
field Kq into constant functions won't necessarily do the same for another Killing field 
K, such that coordinates w.r.t. which every Killing field has constant coefficients (which 
was the content of Theorem 1 in p]) do not necessarily exist. Herewith the reasoning 
leading to the uniqueness of G3 Petrov type I silent universes (Corollary 1 of [?]) breaks 
down. 

Fortunately, there is a related but essentially different way to prove the G3 
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uniqueness result. The key point is that, contrary to gauge dependent quantities such as 
the metric components, invariantly- defined geometric objects must be invariant under 
the isometries f. For example, scalar invariants such as the matter density and the 
eigenvalues of Weyl and expansion tensor must by annihilated by any Killing field, and 
this yields sufficient extra information to prove the key Lemma 3 below. In particular 
for Petrov type I, we note that the eigenvector fields e a of the Weyl tensor and hence the 
corresponding rotation coefficients T abc := e a ■ V c eb are examples of invariantly-defined 
geometric objects and hence are preserved under local isometries. 

2. Mathematical setting 

In what follows, U will always denote an open subset of a spacetime for which the 
conditions in the definition of a silent universe are satisfied everywhere. For such a 
subset, we will say the spacetime is BI if it is an orthogonally spatially homogeneous 
Bianchi I model on U. The set of smooth functions on U will be denoted J~(U), and all 
reasonings will only involve a finite number of polynomial combinations Fi (including 
0) of such functions. Because of this, we can always assume U small enough such that 
for all couples (Fi,Fj) either F»(p) = Fj(p),Wp G U (further denoted by Fi = Fj) or 
Fi(p) + Fj(p)yp G U (denoted by F { ^ Fj). 

For silent universes, the only remaining dynamical variables are the matter 
density /i, the expansion eigenvalues 9 a , the Weyl eigenvalues E a , the purely spatial 
Ricci rotation coefficients 2q a = —e a ■ V a ~i&a-i and 2r a = e a ■ V a+ ie a+ i, and 



m a := d a (n). In these variables, the first contracted Bianchi identity reads cfo/i = —^9 
(with 9 the expansion scalar), and for the derivatives of the expansion and Weyl 
eigenvalues E a , resp. 9 a one has [Oj: 

do9 a = — E a — 9 2 a — — , (2) 
doE a = 29 a+ i x a -i — 29 a _i x a+ \ — ^ a a , (3) 

x a d a 9 a = 6E a (-h a+1 q a + ha-i r a ) - -^m a , (4) 
d a 9 a+1 = -4/i Q _ir Q , d a 9 a -i = -4h a+1 q a , (5) 

d a E a+1 = - 4x a _i r a d a E a ^ = ~4:X a+1 q a (6) 

o o 

Here the expressions have to be read modulo 3 (a running from 1 to 3), o~ a = 9 a — l/3 9 

denote the shear eigenvalues and 2h a := 9 a+1 — 9 a -i, 2x a := E a+ i — E a -\- We 



remark that a particular result stated in these variables actually yields a triple of results, 
by cyclic permutation of the indices. We will use this in Lemma 1(a) and Lemma 3. For 
most of our purposes, it will be favorable to work with the following linear combinations 
of 9 a and E a : 

U 3 TP 3 2 + 03 / 7 x 

xt, hi, ei := —xEi, si := h ■= — ^ — • V) 

f see § 8.4 for an account 
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The expressions for the derivatives of these variables (in terms of themselves and 
q a ,r a ,m a ) can be readily deduced from (J2J)-©- We shall frequently make use of the 
following three results: 

(rl) 9 a = 9p {hi — 0, Si — hx,s\ = —hi) implies E a = E@ {x\ — 0, €\ — Xi,Ei = — Xi), 
see i; 

(r2) if all q a and r a vanish then the spacetime is BI, see [HJ EI] ; 

(r3) the only equilibrium point of the autonomous dynamical system formed by dofi = 
-{i9M and © is (E a = 0,6 a = 0), see 0111]. 

3. Recovering the uniqueness result for G3 silent universes 

Lemma 1 (special cases). If (a) 9 a = (bi = Si,&i = — = hi), (b) bi = or (c) 

^i^i = eihi on U, then the Petrov type is or D. 

Proof, (a) It suffices to prove this for a = 1 (cf. supra). Consecutive time derivatives 
of bi = si, substituting previously obtained equations in each step, yield /i = 4ei, 
— eiSi + xihi = and (ei — a;i)(ei + xi) = 0, i.e. Ei = E 3 or Ei = E 2 . (b) The same 
procedure for bi = yields /i = —Qh\ — 2e\ and Xihi = 0, hence E 2 = E 3 by (rl). (c) 
For xiSi = eihi this procedure immediately yields Xi(si — hi)(si + hi) = 0, so again the 
Petrov type is or D by (rl). □ 

Lemma 2. An algebraically general silent universe U for which si = chi, with c a 
constant function, is necessarily BI. 

Proof. The same procedure as in Lemma 1, applied on si = chi, yields 
ei = (c 2 — l)hf + cxi, hi(c — l)(c + l)(ch\ — bihi + xi) = 0. By (rl), the last equation 
implies x\ = hiipi — chi), and hence ei = hi(cbi — hi) by substition. Two further 
time evolutions yield n = —h\ — 2chibi + 3bf and an identity. Substituting the obtained 
expressions for si, ei, xi, n into ©-(JSJ) yields 3 homogeneous linear systems of 5 equations 
in the variables d a bi,d a hi,q a ,r a ,m a (for each a separately), the coefficients of which 
depend polynomially on bi and hi. The product of the 3 corresponding determinants is 
constantly proportional to h\ 2 (c — l) 8 (c + l) 8 (/ii — &i) 2 (&i + hi) 2 , and this is non-zero 
by (rl) and Lemma 1(a). Hence the linear systems only allow the trivial solution; the 
result follows by (r2). □ 

Lemma 3. When a vectorfield of the form d a + g a do, defined on an algebraically 
general silent universe U (g a G J~(U)) annihilates the scalars \i and (j7j) . then g a = 0. 

Outline of proof. It suffices to prove this for a = 1 (cf. supra), and we set 
X := di + gido- By substituting the annihilation condition mi = <7i /i (3&i — si) for 
/i into X(xi) = X(hi) = X(e-\) = X(xi) = xiX(si) = X(bx) = respectively, one 
gets five equations eqi, . . . , eq§ which form a homogeneous linear system in the variables 
qi,ri,gi, the coefficients of which depend polynomially on /i, x\, hi, e\, si, b\. 

Assume gi ^ 0. This implies deti = det 2 = det 3 = 0, where detj is the 
determimant corresponding to the linear system (egi, eq 2 , eqi), i = 3,4,5. The 
combinations det 2 — g/ixdeti and 6eidet3 — 2sideti can be factorized as (xiSi — eihi)pi, 
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(xiSi — e\hi)p2 respectively, such that p\ = p 2 = by Lemma 1(c). To complete the 
proof we will use three additional relations ^3 = p^ = p§ = 0, where p 3 is defined by 
dopi + (6bx — 2si)pi := 4/iiP3 andp4,ps are the respective polynomials <9oP2 + (5&i — si)f>2 
and d p 3 . 

We scale fi,e\,xi with hi and bi,s± with hi] this comes down to specializing 
hi = 1 in pi, . . . ,p 5 , and we will denote the resulting variables and polynomials with 
a bar. The respective resultants w.r.t. /Z of the couples (Pi,q) with q = p 2 ,p 3 ,p4,p 5 
yield four irreducible polynomial relations U\ = u 2 = 113 = u 4 = for the variables 
£i,xi,bx,si. Next, the resultants w.r.t. ~e~\ of the couples (u^jUi), i = 1,2,4, are 
respectively of the form x\ b\ V\ (xi ,si,b\), x\b l V2 (xi , ~s~i , b\ ) and x\ b l v 3 (afi , s\ , b\) , such 
that V\ — Vi — v 3 = by Lemma 1(b). Finally, the resultants w.r.t. x\ of (^1,^2) and 
(ui,u 3 ) yield 

(&! - sx) 2 ^ + S!) 2 P! p 2 = 0, (Si - s 1 ) 2 g 1 g 2 = (8) 

respectively, where P 1; P 2 , Qi, Q2 and 61 + Si are different and irreducible polynomials 
in bi and s\. By Lemma 1(a) the factor (61 — si) 2 in both left hand sides of (JBJ) can be 
stripped off. Since the remaining expressions have no common factors, their resultant 
w.r.t. bi must yield a non-trivial polynomial relation P(si) = f. Hence ~s~i — j^- is con- 
stant on U, such that the spacetime is BI by Lemma 2. But then di cx -J^j annihilates 
all E a , 9 a and /1, and hence the same is true for gid = (di + gid ) — di, which leads to 
a contradiction with (r3). We conclude that our assumption gi 7^ was false and this 
finishes the proof. □ 

Remark. The purpose of taking the combinations pi,p 2 ,P3 instead of det2,det3 
and dopi from which they are derived, is a lowering of the degree from the start. This 
results in a reduction of the total computation time by a factor of 10 2 — 10 3 . 



Theorem. An algebraically general silent universe U with a group G3 of isometries 
is necessarily a BI spacetime. 

Proof. Since the Petrov type is I, the orthonormal basis {eo, e a } is uniquely 
determined (principal Weyl tetrad); for any point p 6 U, the dimension of the linear 
isotropy group of p is and hence the dimension of the orbit O p under the supposed 
G3 isometry group is 3 §. This implies the existence of three functionally independent 
Killing vectorfields Ki on U, i = 1 . . . 3. Writing the expansions Ki = Kf do + Kf d a = 
9 a w.r.t. the principal Weyl tetrad in block matrix form: 

do 



K 



Ki 
K 2 
K 3 



K°d + A 



0, 



M 



di 

2 

0, 



(9) 



\ see e.g. ^21, PP- 159- More explicitly, this can be checked by specializing e.g. s% — 2 in the remaining 
expressions: the degrees of b\ don't drop, and calculating the specialized resultant gives a non-zero 
number, such that the original resultant cannot be identically (see also ^0] for this argument). 
§ see CP, § 8.3, 8.4 and 9.2 for an account 
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where K° G JF(£f) 3xl , A G T{Uf xZ and M G ^(tf) 3 * 4 , we thus have rank(M)= 3 and 
hence rank(v4)= 2 or 3 everywhere on U. 

(i) rank(A)=2. 

In this case, there exist functions fi, f 2 , f3 G ^F(U) such that Y^i=i fi ^i* = 0, where A itf 
denotes the i th row of A. Applying the same linear combination to (|SJ) we get 

3 3 
i=l »=1 

where the coefficient tp nowhere vanishes on U because rank(M)=3. Hence, since the 
Ki annihilate all Weyl and expansion eigenvalues, the same is true for do, which leads 
again to a contradiction with (r3). Hence this case cannot occur. 

(ii) rank(A) = 3. 

In this case A is everywhere invertible on U. When applying the inverse, say C, on the 
left of (jOJ, the a-component of the resulting equations yields 

d a + g a d = C ai Ki (10) 

with g a := C a iK^. Hence d a + g a do (a = 1,2,3) annihilate all scalar invariants, and 
this implies g a = by Lemma 3. Thus all d a annihilate (in particular) 9 a+ \ and 9 a -i, 
which forces all r a and q a to vanish by (jSJ) and (rl); the result follows by (2). □ 
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